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Abstract. Charged particle multiplicities in A+A and p(p) + p collisions as a function of pseudorapidity, 
centrality and energy are studied in both the nucleon and the constituent quark frameworks. In the present 
work, the calculation of the number of nucleon and constituent quark participants using the nuclear overlap 
model takes into account the fact that for the peripheral A+A and p+p collisions can not be smaller 
than two. A striking agreement is seen between the particle density in A+A and p(p) + p collisions, 
both at mid-rapidity and in the fragmentation regions, when normalized to the number of participating 
constituent quarks. The observations presented in this paper imply that the number of constituent quark 
pairs participating in the collision controls the particle production. One may therefore conjecture that 
the initial states in A+A and p(p) + p collisions are similar when the partonic considerations are used in 
normalization. 

PACS. 25.75.-q, 25.75.Dw, 25.75.Nq, 24.10.Jv 



1 -Introduction 

Quantum Chromo-Dynamics (QCD) is the theory of strong 
interactions. From a theorist's point of view strong inter- 
action means that all dimensionless quantities of the the- 
ory are of the order of unity. Meanwhile, there are at least 
two fundamental facts about strong interactions showing 
that actually it is not exactly so: there are still certain 
small parameters hidden. One is that nuclei can be viewed 
as made of weakly bound nucleons - in a true strong in- 
teractions case they would rather looke like a quark soup. 
The second puzzle is that the nucleon itself can be viewed 
as built of three constituent quarks (objects quasi-free non 
relativistic) - in a true strong interactions case it would 
rather look like a pack of an indefinite number of quarks, 
antiquarks and gluons. The presence of just three discrete 
objects in a nucleon is non-trivial from the point of view 
of the parton model, which permits any large number of 
point-like objects, i.e. quark-partons, inside a fast moving 
nucleon. The two notions can be reconciled if one accepts 
a nucleon with two radii: a nucleon radius R n that deter- 
mines a mean distance between the constituent quarks, 
and a proper radius of the constituent, r q . In such a pat- 
tern, a fast moving nucleon is a viewed as a system of three 
clouds of partons, each containing a valence quark, a sea of 
the quark-antiquark pairs, and gluons. QCD calculations 
support the statment that inside a nucleon there are three 
objects of size of 0.1-0.3 fm £Q. Similar idea is applied in 
the study of multihadron production processes in different 
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types of collisions in the framework of the picture based 
on dissipating energy of participants [2] ■ 

In this paper, I illustrate a modified version on how 
to calculate the number of nucleon and constituent quark 
participants. One problem in previous calculations 3,4 
was confusion in applying the appropriate minimum-bias 
(Npart) to the p(p) + p collisions. A second, more general 
problem in the previous calculation which can be found 
in the Refs. 013115] stems from the fact that in the nu- 
clear overlap model in peripheral collisions the number of 
nucleon and constituent quark participants ((Np art ) and 
(Npart) ) can De smaller than two - this is because, in pe- 
ripheral limit, the overlap integral has a meaning of 1/2 
times the probability to have (Np art ) =2. In the previous 
calculation 3,4,5J, the problem can be found in the es- 
timate of (Np art ) for the peripheral A+A collisions and 
easily from the calculation of the number of constituent 
quark participants for p+p collisions such plotting (Np art ) 
versus collision centrality or impact parameter. A similar 
problem with the calculation of (Np art ) for d+Au colli- 
sions using an optical approch has been reported by D. 
Kharzeev et al. jEj- The goal of the present paper is to 
study the charged particle multiplicities in the nucleus- 
nucleus (A+A) and nucleon-nucleon (p(p) + p) collisions 
in both the nucleon and the constituent quark frameworks. 

The motivations of this study are; nuclear collisions 
are studied at unprecedented energies at the Relativistic 
Heavy Ion Collider (RHIC), revealing several new phe- 
nomena embedded in a large amount of high-quality data 
[7|. However, several aspects of the results on charged par- 
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tide multiplicities in the nucleon framework (results scaled 
to the number of nucleon participants) are still not well 
understood such as; 

1. the charged particle densities near mid-rapidity region 
in Au+Au at 200 GeV depend strongly on the colli- 
sion centrality. The densities have been plotted as per 
nucleon participants pair 8 in order to disentangle 
pure nuclear effects. The common explanation of the 
phenomena is due to the hard processes. 

2. the charged particle density at mid-rapidity region from 
A+A collisions are substantially higher than those of 
p(p) + p collisions at the same energy 0E3], 

3. the charged particle densities at the fragmentation re- 
gion from A+A collisions are substantially higher than 
those of p(p) + p collisions at the same energy 0], 

4. the integrated total charged particle in A+A collisions 
as a function of number nucleon participants is higher 
than p(p) + p collisions at the same energy indicating 
that there is no smooth transition between A+A and 
nucleon- nucleon collisions 11 . 

The investigations on the aspects 1. and 2. have been al- 
ready started in Refs. P]|4*||^] using the old version which 
contains problem in the calculations of (Np art ) and (Np art ) 
for peripheral A+A collisions and p(p) + p collisions. In 
this paper, I will review the aspects 1. and 2. and I will 
extend the investigation to study the aspects 3. and 4. in 
both the nucleon and constituent quark frameworks, using 
the new version of the calculation. 

2 -Calculation of the Number of Participants 

The number of nucleon participants noted by (Np art ) , 
and the number of constituent quark participants noted 
by (Np art ) , are estimated using the nuclear overlap model 
in a manner similar to that used in Refs . 3,4,5 but I intro- 
duce in the present work a modification of the calculation 
procedure taking into account that for the peripheral col- 
lisions A+A and for p+p collisions in both framewroks; 
the (Np art ) and (Np art ) can not be smaller than two. The 
nuclear density profile is thus assumed to have a Woods- 
Saxon form, 

Mr) = 1 , l u° 5 w, i ; (!) 



1 



exp[(r - R n )/d] 



where no is the normal nuclear density, R n is the nucleus 
radius and d is a diffuseness parameter. 

For nucleus- nucleus (A+B) collisions, the number of 
nucleon participants, (Np art ) , is in the present work cal- 
culated using the relation, 
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Fig. 1. Quantitative evaluation of the model calculations ex- 
pressed as the ratio of the average number of nucleon partici- 
pants of the PHOBOS Glauber calculations to the present 
work as a function of centrality. Panel a), b) and c) correspond 
to Au+Au collisions at ^/s~^= 19.6, 62.4 and 200 GeV, re- 
spectively. The gray bands correspond to the systematic errors 
on the (Np art ) of PHOBOS Glauber calculations. 
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where T(b) = J_°° dzn A {\/b 2 + z 2 ) is the thickness func- 
tion and P AB {b) is defined by: 



Fig. 2. Ratio of (Np art ) /(Np art ) obtained from the present 
work presented as a function of the number of nucleon par- 
ticipants for Au+Au collisions at ^/s NN = 19.6, 62.4 and 200 
GeV. 
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Fig. 3. Distributions of constituent quark participants, (Np art ) 
, obtained from the present work for p+p at v /s JVJV = 53, 200, 



900 GeV represented as a function of collisions centrality. 



Table 1. Numbers of constituent quark participants, (Np art ) 
, obtained from the present work elucidated as a function of 



collision centrality in Au+Au collisions at .Vs^ 
and 200 GeV. 



7= 19.6, 62.4 



Uentranty \ /o) 






i n ^\ T 

ly.o ueV 


Bin 








0-3 


952.5 


907.2 


869.4 


3-6 


837.7 


796.0 


761.4 


6-10 


731.4 


693.1 


661.8 


10-15 


614.5 


580.4 


552.8 


15-20 


501.7 


472.1 


448.4 


20-25 


408.5 


382.9 


362.7 


25-30 


331.6 


309.7 


292.6 


30-35 


265.3 


246.8 


232.4 


35-40 


209.1 


193.5 


181.6 


40-45 


161.6 


148.8 


139.1 


45-50 


123.1 


112.8 


105.1 



where Tab (6) is the overlap function defined as the prod- 
uct of the thickness functions of the colliding nuclei A and 
B, integrated over the two transverse dimensions: 

T AB (b) = J <? S 'T A (s')T B {s' - b) (4) 

A and B are the mass number of the two colliding 
nuclei and the <J™fj is the inelastic nucleon-nucleon cross 
section. In the present work I use the FRITIOF param- 
eterization R n = 1.16 A 1 ' 3 - (1.16) 2 A- 1 / 3 fm and d = 

0. 54 fm which are well within the measurements of elec- 
tron scattering from Au nuclei |12| . 

Fig. OD shows the ratio of (N£ art ) obtained by PHO- 
BOS ^31 an d to the (Np art ) from the present work for 
Au+Au at y/s^= 19.6, 62.4 and 200 GeV. The (N£ art ) 
and (Np art ) from the present work are presented on Ta- 
ble for Au+Au collisions at RHIC energies. 

The number of constituent quark participants, (Np art ) 
, is calculated in a similar manner by taking into account 
the following changes related to the physical realities and 
also taking into account that in peripheral collisions the 
number of (Np art ) can not be smaller than two: 

1. the density is three times that of nucleon density with 
n q = 3n = 0.51 fm" 3 , 

2. the cross sections <jcq — er™^/9, 

3. the mass numbers of the colliding nuclei are three times 
their values, keeping the size of the nuclei same as in 
the case of (N£ art > . 

Fig. Elpresents the ratio of (Np art ) /(Np art ) as a func- 
tion of (Np art ) for Au+Au collisions at RHIC energies. 
The ratio shows that the correlation between (Np art ) and 
(Np art ) is not linear and that it depends on the colliding 
energy. 

For proton-proton (p+p) collisions the same procedure 
has been used to calculate the number of constituent quark 
participants by using A = 3 and B = 3 and the nuclear 
density profile is assumed to have a sharp sphere form 
with uniform radii of 0.8 fm [T3|. The reason to use sharp 
sphere density profile is because the Wood-Saxon density 



Table 2. Values of overlap function TAB and TA for impact 
parameter b=0 calculated using densities profiles Wood-Saxon 
and sharp sphere for p(A=l)+B collisions using the nuclear 
overlap model 



profile 


A 


B 


TAB(0)/TA(0) 


Wood-Saxon 


1 


236 


0.9 


Wood-Saxon 


1 


12 


0.6 


sharp sphere 


1 


236 


1.0 


sharp sphere 


1 


12 


0.9 



Table 3. Values of (Np art ) obtained from the present work for 
most central and minimum-bias (min— bias) of p(p) + p colli- 
sions at several collisions energies. 

ysTJT: (GeV) I 53 I 200 I 540 I 630 I 900 I 1800 

ocq = (Jiviv/9 3.89 4.66 5.33 5.44 5.66 6.22 

(Ng art ) (0-6%) 4.37 4.94 5.46 5.55 5.72 5.98 

(Ng art ) (min-bias) 3.10 3.29 3.47 3.51 3.57 3.73 



profile becomes unrealistic for low A. The systematic error 
related to this approximation can be estimated by run- 
ning the code with B=l (for p+A), and comparing ob- 
tained TAB(0) value with TA(0) (for impact parameter 
b=0). The result is presented in Tabled The discrepancy 
is larger for Wood-Saxon profile with low A. The (Np art ) 
for central, 0-6%, and minimum bias of p+p collisions are 
presented on the Table 01 Fig. |3 shows the distributions 
of the number of constituent quark participants obtained 
from the present work elucidated as a function of collisions 
centrality for nucleon-nucleon (p+p) collisions at y / s NN = 
53, 200 and 900 GeV. 



3 -Physics Results and Discussions 

Fig.EJshows the centrality dependence of the charged par- 
ticle density per participant pair. Fig.^Ja shows the PHO- 
BOS results PT3] on dN ch /dr]\ lv <! per nucleon partici- 
pant pair for Au+Au at 19.6, 62.4 and 200 GeV. The 
centrality dependence of the mid-rapidity yields has often 
been interpreted in a two component picture of particle 
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Fig. 4. Particle density, dJVch/d*7||n|<ij divided by a) the num- 
ber of nucleon participants pair, b) the number of constituent 
quark participants pair in Au+Au and p(p) + p collisions at 
three energies [8115120). The dashed lines in panel b) corre- 
spond to the average of Au+Au data at given energy. The 
p(p) + p for non single diffractive (NSD) data are represented 
with similar symbols as Au+Au collisions at the same energy 
at (Np art ) —2. The star symbols in panel a) and b) correspond 
to p(p) + p inelastic data at ^/s NN = 200 GeV. The p(p) + p 
data in in panel b) have been normalized by the minimum-bias 
(N^ art ) . 



production by soft and hard processes. As the beam en- 
ergy increases, particle production from hard processes, 
which exceed the number of participants pairs by a fac- 
tor ~ 5-6 in central events for y/s NN ranging from 19.6 to 
200 GeV, is expected to dominate over that from soft pro- 
cesses as the mini-jets cross sections increase |15|. Fig.QJb 
shows the centrality dependence of dNch/drjl^^i per con- 
stituent quark participants pair. I observe a constant or a 
slightly decreasing dependence of (dN/drj) /((Np art ) /2) on 
centrality from peripheral to central collisions for Au+Au 
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Fig. 5. Particle density per constituent quark participant pair 
(open symbols) and particle density per nucleon participant 
pair (solid points) produced in central (0-6%) nucleus-nucleus 
(A+A) collisions presented as a function of collision energy at 
AGS, SPS 1161171181 and RHIC [T§] and in p(p) + p collisions 
(inelastic) [U]. The errors bars correspond to the systematic er- 
rors. The solid line represents a linear fit through solid symbols 
for A+A data, f A A = -0.287 + 0.757rn(Vs)- The dashed dot- 
ted line corresponds to the fit through solid symbols of p(p) + p 
collisions, f pp = 2.25 - 0.411n(^s) + 0.091n 2 (^s). The dashed 
line corresponds to a linear fit trough the open symbols in the 
constituent quarks framework, f p (p) p /AA — —0.06 + 0.31n(y / s). 
It should be noted that the dN c h/dr]\ i^iof p(p) + p data 
(open symbols) have been normalized by minimum-bias (Np art ) 
presented in tabled 



at 200 GeV and good scaling of {dN / drf) / {{^ axt ) /2) on 
centrality for Au+Au at 19.6 and 62.4 GeV. In Fig.Hb, I 
extend the study presented in Ref. [Hj from Au+Au at 200 
GeV to 19.6 and 62.4 GeV using the new calculations. I 
agree with the interpretation presented in Ref. that the 
experimentally observed increase of (dN^/dr)]^^)/ ((Np art ) 
/2) can be explained by the relative increase in the number 
of interacting constituent quarks in more central collisions. 
In the present work, I add to this study dNch/drj^^x of 
p(p) + p inelastic as well NSD collisions. In Fig.|Ub, I ob- 
serve that the p(p) + p data are in good agreement with 
Au+Au data for different collision centrality at the same 
energy. It should be noted that in Fig.^b the dN c h/ d-q^^i 
of p(p) + p data (open symbols) have been normalized by 
minimum-bias (Np art ) presented in table 01 
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Fig. 6. Integrated total charged particles obtained from 
Au+Au collisions and p(p) + p inelastic as well NSD collisions 
at the same energy ys^7= 200 GeV (8113120) normalized by a) 
number of nucleon participant pair, b) number of constituent 
quarks participant pair. The errors shown with vertical bars 
are full systematic errors. In panel b) the N c h for p(p) + p 
data (open symbols) have been normalized by minimum-bias 
(Np art ) . The dashed lines are guide lines relatively to the most 
central collisions in Au+Au system 



Fig. shows the primary charged particle density for 
central collisions at mid-rapidity divided by a) the number 
of participant nucleon pairs ((Np art ) /2) as solid symbols 
and b) the number of participant constituent quark pairs 
((Npart) /2) as open symbols presented as a function of col- 
lision energy. The data are for Au+Au collisions at AGS, 
Pb+Pb collisions at the CERN-SPS [HI17lllcl| and for 
Au+Au and Cu+Cu collisions at RHIC |19) . Also shown 
for comparison are results from p(p)+p collisions |2()| . 

In the nucleon participants framework (solid symbols 
in FiglSJ), the particle density per nucleon participant pair 
for A+A collisions (solid points) shows an approximately 
logarithmic rise with ^/s^ over the full range of collision 
energies. The comparison of the particle density per nu- 
cleon of Au+Au to Cu+Cu collisions at the same energies, 
^/snn = 62.4 and 200 GeV indicates that in symmetric 
nucleus-nucleus collisions the density per nucleon partici- 
pant does not depend on the size of the two colliding nuclei 
but only on the collision energy. This means that for Si+Si 
collisions at ^/snn = 200 GeV, the particle density per 
nucleon participant will be similar to Au+Au collisions at 
the same energy. Also I observe that the charged particle 



multiplicity per participant nucleon pair (solid symbols) 
in A+A collisions is higher compared to p(p) + p colli- 
sions at the same energy. It should be noted that number 
of nucleon participants for p(p) + p has been chosen to be 
2 «N£ art ) ( PP ) =2). 

In contrast, I observe that the particle density per con- 
stituent quark participant pair (open symbols in FigEJ is 
similar for nucleus-nucleus collisions and nucleon-nucleon 
collisions at the same energy. It thus appears that using 
partonic participants accounts for the observed multiplic- 
ity in both A+A and p(p) + p collisions. One may there- 
fore conjecture that the initial states in A+A and p(p) + p 
collisions are similar. It should be noted that Figgis sim- 
ilar to Fig. 6 presented in Ref. : 3. but it was done indepen- 
dently, extended to the Cu+Cu system, and used the new 
version of calculation presented in this paper. 

Fig. shows the integrated total charged particle di- 
vided by a) per nucleon participant pair (A r c ^/((Np art ) 
/2)), b) per constituent quarks participant pair (N c iJ ((Np art ) 
/2)) for Au+Au and p(p) + p collisions at the same en- 
ergy, 200 GeV. In the nucleon participant framework, panel 
a), the integrated total charged particle in Au+Au col- 
lisions, as a function of number nucleon participants, is 
higher than p(p)+p collisions indicating that there is no 
smooth transition between A+A and nucleon-nucleon col- 
lisions This observation can be strongly confirmed 
by adding the preliminary data from Cu+Cu presented 
by PHOBOS collaboration JHj (not shown in this figure) 
which correspond to the peripheral data of Au+Au system 
at the same energy. These comparisons between A+A and 
p(p) + p collisions are, however, based on scaling with the 
number of nucleon participants. In the constituent quarks 
framework, panel b), the A^ c fc/((N^ art ) /2) of p+p NSD 
(or inelastic) collisions agree with the central Au+Au col- 
lisions at the same energy. 

In general, the charged particle production in the frag- 
mentation region is thought to be distinct from that at 
mid-rapidity, although there is no obvious evidence for 
two separate regions at any of the RHIC energies. This 
observation is made based on the dN c h/drj distributions of 
charged particle presented in Ref. [TUll2*T] . Fig.[7Jshows the 
charged particles produced in the fragmentation region for 
the most central (0-6%) Au+Au collisions at four RHIC 
energies compared to p(p) + p (inelastic and (NSD)) colli- 
sions at 200 GeV. When normalized to (N^ art ) /2, Fig.|7Ja, 
I observe that the multiplicity in the limiting fragmenta- 
tion region in A+A collisions is higher than for p(p)+p 
collisions at the same energy, y/s NN =200 GeV. If, however, 
the comparison is carried out for multiplicities normalized 
to (N^) /2, FigE|b, A+A and p(p) + p collisions exhibit 
a striking degree of agreement. Again, this observation im- 
plies that the number of constituent quark pairs partici- 
pating in the collision controls the particle production in 
the central collisions. 



4 -Conclusions 

The charged particle production results from A+A and 
p(p) + p collisions have been compared based on the num- 
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Ml " Ybeam 

Fig. 7. Pseudorapidity distributions of charged particle for 
Au+Au collisions at RHIC energies |5| and p(p)+p collisions at 
200 GeV |20|. The distributions have been shifted to r\ — ybeam 
in order to study the fragmentation regions in one of the nu- 
cleus rest frame. Panels a) and b) correspond to dN^/drj distri- 
butions scaled to the number on nucleon participants pair and 
to the number of constituent quark participants pair, respec- 
tively. The systematic errors have been shown just for Au+Au 
at 200 GeV for clarity. In panel b) the p(p) + p data have been 
normalized by the minimum-bias (Np art ) . 



ber of nucleon participants and the number of constituent 
quark participants. In both normalizations, I observe that 
the charged particle densities in Au+Au and Cu+Cu col- 
lisions are similar for both ^/s JI „= 62.4 and 200 GeV. 
This implies that in symmetric nucleus-nucleus collisions 
the charged particle density does not depend on the size of 
the two colliding nuclei but only on the collision energy. In 
the nucleon participants framework, the particle density 
at mid-rapidity as well as in the fragmentation region from 
A+A collisions are higher than those of p(p) + p collisions 
at the same energy. Also the multiplicity of total charged 



particle, in A+A collisions, as a function of number nu- 
cleon participants is higher than p(p) + p collisions at the 
same energy indicating that there is no smooth transition 
between peripheral A+A and nucleon- nucleon collisions. 
However, when the comparison is made in the constituent 
quarks framework, A+A and p(p) + p collisions exhibit a 
striking degree of agreement. The observations presented 
in this paper imply that the number of constituent quark 
pairs participating in the collision controls the particle 
production. One may therefore conjecture that the initial 
states A+A and p(p) + p collisions are similar when the 
partonic considerations are used in normalization. 
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